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The mode number of the Dirac operator scales with an exponent related to the mass anomalous 
dimension 7,,,. This relation holds both in IR-conformal systems, as well as in confining systems 
for large enough eigenvalues. We investigate the Nf ~ 4, di and 12 flavor SU(3) systems at sev- 
eral couplings near the chiral limit, and show that in general the scaling exponent varies with the 
eigenvalue, describing the dependence of Ym on the energy (or, equivalently, on the running cou- 
pling). This energy dependence can be explored even with fixed lattice parameters (bare coupling 
and mass). We find that for the 4 flavor system the mass anomalous dimension decreases as the 
energy increases, consistent with perturbative expectations. For the 8 flavor system the energy de- 
pendence is too weak to be observable at present. The 12 flavor system at our strongest couplings 
shows the anomalous dimension increasing with energy, consistent with backward flow and the 
presence of an infrared fixed point. At weaker couplings we determine a preliminary value for the 
mass anomalous dimension of the 12 flavor system at the infrared fixed point, ■j^, = 0.27(3). 



X 
S3 



30th International Symposium on Lattice Field Theory 
24-29 June 2012 
Cairns, Australia 



* Speaker. 



© Copyright owned by the author{s) under the terms of the Creative Coinmons Attribution-NonCoinmercial-ShureAHke Licence. 



http://pos.sissa.it/ 



Mass anomalous dimension from Dime eigenmode sealing 



Anna Hasenfratz 



1. Introduction 

The lattice study of gauge systems with many flavors or higher representation fermions is 
motivated in part by the possibility of electroweak symmetry breaking via new strong dynamics. 
The discovery of a 125GeV Higgs-like particle puts severe restrictions on many BSM scenarios, 
including models of new strong dynamics, which are now required to predict a light composite 
Higgs. Even if the particle discovered at 125GeV turns out to be an elementary scalar, the study of 
these strongly coupled systems remains important not only for theory but model building as well. 

While the numerical simulation of many-fermion systems is not fundamentally different from 
well understood QCD simulations, it is increasingly clear that QCD-like analyses are not always 
optimal. The large number of fermions requires working with relatively strong couplings where 
unusual phases, unexpected lattice artifacts and spurious UV fixed points can influence the critical 
behavior [1]. We use several complementary techniques and approaches to investigate these lattice 
models, comparing and contrasting systems with different numbers of fermions to distinguish con- 
formal and chirally broken, confining behavior. In this paper we consider the scaling of the Dirac 
eigenmodes; our work on the phase structure, finite temperature behavior, and Monte Carlo renor- 
malization group (MCRG) analysis is discussed in other contributions to these proceedings [2, 3]. 

Our gauge action consists of fundamental and adjoint plaquette terms and we use nHYP 
smeared staggered fermions. The details of the lattice action can be found in Ref. [4]. In this 
work we are interested in the weak coupling behavior of these systems, and we avoid the single- 
site shift symmetry broken (S^) lattice phase discussed in Refs. [1, 2]. The Nf = 4, 8, 12 and 
16 flavor systems have been investigated extensively by several other groups; recent references 
include [5, 6, 7, 8] and earlier works are reviewed in Ref. [9]. 

2. The eigenmode density and mode number 

The eigenmodes of the Dirac operator contain a wealth of information about the dynamics 
of lattice systems. When the infrared behavior is captured in a Random Matrix Theory (RMT) 
universality class, the distribution of the low eigenmodes predicts physical quantities like the chiral 
condensate Z. Conformal theories are more difficult to analyze as there are no RMT predictions for 
the individual eigenmodes. One approach is analyze the low-lying eigenmodes on several different 
volumes; a simple finite-volume scaling fit predicts the mass anomalous dimension, though with 
large systematic errors [10, 1]. 

It is more reliable to consider the eigenmode density p(A) or its integral, the mode number 



where V is the volume of the system. The mode number is renormalization group invariant even 
at finite fermion mass, and can be used to predict the chiral condensate in QCD-like systems with 
high precision [11]. 

In conformal systems where the low energy dynamics is governed by an infrared fixed point 
(IRFP), the general scaling form p(A) oc A" predicts 
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where L = V^^^ is the linear system size. Using the renormalization group invariance of the mode 
number, this relates the exponent a and the scaling dimension of the mass ym 

4 

3'm = l + 7«. = — j-- (2-3) 

Ref. [12] used the mode number, evaluated with a stochastic method [11] in a wide eigenvalue 
range, to predict the mass anomalous dimension for SU(2) gauge theory with two adjoint fermions. 
In this work we follow a similar approach both in chirally broken and IR-conformal theories. 

3. The energy dependence of the scaUng dimension 

The eigenvalues A of the Dirac operator have dimension of mass. From A ~ to the cutoff 
scale A ~ a^^ the eigenvalue distribution reflects the properties of the system from the infrared to 
the ultraviolet. In a chirally broken system we can distinguish three regions: 

I. Low energy region, below the chiral symmetry breaking scale. Here the eigenvalues describe 
the infrared behavior of the system, as is reflected in the Banks-Casher formula p{0) = 2L/7l 
[13]. The correction in A can be calculated in chiral perturbation theory and the scaling of the 
mode number predicts the chiral condensate £. 

II. Above the chiral symmetry breaking scale but still well below the cutoff the system is gov- 
erned by the Gaussian fixed point and the eigenvalue density scales according to Eq. 2.2 with 
an energy dependent scaling exponent. At weak coupling, perturbation theory predicts the 
one-loop universal running of the anomalous dimension as Y,„ = 6Cfg^/{47l)^ + C?(g^) for 
fermions in the fundamental representation. 

III. The high energy, ultraviolet part of the spectrum for lattice systems is non-universal, domi- 
nated by discretization effects. 

In chirally broken systems the anomalous dimension 7,,, depends on the running coupling and 
consequently the energy scale. This is reflected in the scaling form of Eq. 2.2 as a scaling exponent 
that depends on A as a measure of the energy scale. Just above the chiral symmetry breaking scale 
7,,, is large; in the standard folklore 7,, 1, therefore a ^ I. At high energies perturbation theory 
predicts 7m \ and a /^3. 

Conformal systems have only regions II and III. The low energy end of the spectrum reflects 
the properties of the infrared fixed point and the scaling form Eq. 2.2 predicts the anomalous di- 
mension at the IRFP, 7,, = 7];,. Starting from a weak bare coupling, 7,, at high energies will again 
approach the perturbative 7,, = 0. There is no perturbative prediction for what happens if the bare 
coupling is in the strong coupling side of the basin of attraction of the IRFR ^ 

The above discussion might sound speculative at this point. In the following we will support 
it using data from our 4, 8 and 12 flavor simulations. 



The location of the IRFP depends on the specific renormalization group transformation considered; it is not a 
physical observable, and the strong coupling side of the IRFP is not well-defined. However it is possible to observe 
whether or not a lattice system is in the strong coupling side of the IRFP's basin of attraction. 
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Figure 1: The eigenvalue density p (A ) at j3/r = 2.8 of the 12 flavor system. Left panel: 16^ x 32 volumes 
at various sea fermion mass values. Right panel: mass am — 0.0025 on various volumes. 

4. Finite mass and finite volume effects 

The scaling form of Eq. 2.2 assumes infinite volume and vanishing fermion mass. In practice 
neither of these conditions are satisfied in lattice simulations. The authors of Refs. [11, 12] argue 
that the fermion mass affects only the low eigenmodes. Since we use staggered fermions we eval- 
uate the eigenvalues of the massless Dirac operator, on configurations generated with nonzero sea 
fermion masses. The left panel of Fig. 1 shows p(A) for different sea fermion masses in the 12 
flavor system at j3f = 2.8 on 16^ x32 volumes. This coupling is safely on the weak coupling side of 
the^ph ase [1]. For am > 0.02, the eigenvalue density depends strongly on the mass, and appears 
unlikely to become mass independent even at larger A. On the other hand for am < 0.01 the mass 
dependence rapidly disappears as A increases, suggesting that here it is possible to reach the chiral 
Umit by simple extrapolation. 

The right panel of Fig. 1 illustrates the volume dependence of p (A ) for the smallest mass in 
the left panel, am = 0.0025. The system is chirally symmetric on all four volumes considered, 
with p(0) = 0. At the small A where the density becomes nonzero, there is transient volume 
dependence. All four volumes produce consistent results for larger A that are still well below 
the cutoff scale. While there remains a small volume dependence even in this regime, an infinite 
volume extrapolation is feasible. 

Even in systems known to exhibit spontaneous chiral symmetry breaking we find consistent, 
largely volume independent p(A) for A above the chiral symmetry breaking scale. By comparing 
data on different volumes and at different sea fermion masses it is possible to identify the regime 
where both the infinite volume and chiral extrapolations are feasible. 

5. Scaling of the mode number 

In this section we present results for the mass anomalous dimension based on the scaling 
of the mode number (Eq. 2.2) within region II. We use 12^ x24, 16^ x32 and 24^ x48 volumes. 
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Figure 2: The anomalous dimension 7,,, based on the scaling relation of Eq. 2.2 as a function of X, on 
volumes 12^ x24 (red squares), 16^* x 32 (green circles) and 24^ x48 (blue triangles). Left panel: Nf ^ A 
flavor system at j3f = 6.6, am = 0.0025. Right panel: Nf ~% flavor system at /3f = 4.8, am = 0.0025. 



generate ensembles of 50-100 thermalized configurations, and calculate 300-1000 eigenvalues on 
each configuration. The data we present are obtained with mass am = 0.0025, and in some cases 
we consider smaller values to check for finite mass effects. Since we would like to see if and how 
the scaling exponent changes with A, we fit the mode number over a range AA, as 

v(A + AA /2) - v(A - AA /I) = c [(A + AA/2)«+' - (A - AA /2)«+'] . (5.1) 

We carry out this two-parameter fit for fixed volume, mass and gauge coupling; it is fairly stable if 
AA is not too small. 

The left panel of Fig. 2 shows the anomalous dimension = [4/(a + 1) — 1] as a function 
of A in the 4 flavor system at j3f = 6.6. At this coupling and mass all three volumes are chirally 
symmetric, p(0) = 0. The transient finite volume effects at small A are clear, but as A increases 
the predictions from all volumes fall onto the same curve. The interesting feature of this plot is the 
steady decrease of with increasing A, i.e. towards the ultraviolet. This is the qualitative energy 
dependence expected from the perturbative relation = 6Cfg^/(47r)^. With g^(A) determined 
from the one-loop beta function, this relation gives a good overall description of the data, though 
with coefficients that differ from their one-loop values by a factor of three. At small A , 7,,, ~ 1 , and 
indeed at this coupling the 24^* x 48 system appears almost chirally broken. 

The right panel of Fig. 2 shows the A dependence of the anomalous dimension for the 8 flavor 
system at j8f = 4.8. At am = 0.0025 all three investigated volumes are chirally symmetric, with 
p (0) = 0. Unlike the Nf = 4 case, the common 7,,, at which all three volumes eventually settle 
shows very little dependence on A. As expected, the 8-flavor anomalous dimension runs much 
more slowly than for Nf = 4. We find similar behavior at other bare couplings, though the value 
of Yin decreases steadily as /3f moves to weaker coupling. We are currently calculating the running 
coupling in this system, which will allow us to combine results at several jSf into a function Ym{g^) 
that we can compare to perturbation theory. At this point we can only conclude that the 8 flavor 
system exhibits slow running compared to Nf = 4, with large 7,„ ~ 1 at j3/r = 4.8; we cannot yet 
quantify the energy range where this behavior persists. 
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Figure 3: The anomalous dimension Ym based on the scaUng relation of Eq. 2.2 as a function of X for 
the Nf ~ 12 flavor system on volumes 12^ x24 (red squares), 16^ x32 (green circles) and 24^ x48 (blue 
triangles). Left panel: j5f = 3.0, am 0.0025. Right panel: )3f = 6.0, am = 0.0025. 



Fig. 3 shows two representative samples from our Nf = 12 investigations. The left panel is the 
anomalous dimension at j3f = 3.0, a fairly strong coupling but still on the weak coupling side of 
the ^phase. Even at this fairly strong coupling, all three volumes are chirally symmetric, and they 
predict a consistent Ym at A large enough for finite volume effects to die out. However, the energy 
dependence is the opposite of what we observed with Nf = 4: y,,, increases with increasing A. This 
is not consistent with the perturbative behavior at the = Gaussian fixed point, and we interpret 
this result as indicating that the coupling j3f = 3.0 is in the strong coupling side of the basin of 
attraction of the conformal infrared fixed point. 

At weaker bare coupling (larger fif) we observe less energy dependence, as the right panel 
of Fig. 3 illustrates at j3f = 6.0. The finite volume effects are more significant at this coupling, 
but the three volumes settle at a consistent value « 0.3 that is largely A -independent. We find 
similar behavior at other Pf > 6. Considering 6 < jSf < 8, we determine the preliminary result 
= 0.27(3) for the mass anomalous dimension at the IR fixed point. We hope to obtain a more 
precise result for y^^^ in the future by performing a global fit over all the volumes, j3f and ranges of 
A that produce consistent predictions. 

While our result = 0.27 (3) is significantly smaller than that obtained from finite-size scaling 
by Ref. [6], it is compai^able to some of the results reported by Ref. [5], which found different 0.2 < 
)^ < 0.4 depending on the observable used in the finite-size scaling analysis. We also reported a 
larger value from investigations at stronger couplings [1]. The fact that we need jSf > 6 to obtain 
scale-independent results could explain some of the inconsistencies at stronger couplings. 



6. Conclusion 

We presented preliminary results from our analysis of the mode number of the Nf = 4, 8 and 
12 flavor SU(3) models. We argued that a simple scaling form that depends on the mass anomalous 
dimension describes the mode number in an intermediate spectral range even for chirally broken 
systems. Our data show the running of /,„ with the energy in the 4 flavor system, while with 8 
flavors 7„, does not change significantly over the accessible range of scales. Our 12-flavor data are 
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consistent with the existence of an infrared fixed point and show backward running at strong bare 
coupling. At weaker couplings we obtain a preliminary prediction y^^^ = 0.27(3) for the Nf = 12 
mass anomalous dimension at the IR fixed point. 

The Dirac eigenmodes offer an alternative method to study the infrared dynamics of strongly 
coupled systems. Our analysis at present is limited by the number of eigenvalues we can calculate 
on the larger volumes. The stochastic estimator proposed in Ref. [11] and used by Ref. [12] may 
greatly increase the reach of this approach. 
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